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Context. A population of unidentified gamma-ray sources is forming a structure resembling a halo around the Galactic 
center. These sources are highly variable, and hence they should be associated with compact objects. Microquasars are 
objects undergoing accretion with relativistic jets; if such an object has a low-mass, evolved, donor star, it might be 
found in the Galactic halo. If these low-mass microquasars can generate detectable gamma-ray emission, then they are 
natural candidates to account for the halo high-energy sources. 
Aims. We aim to construct models for high-energy emission of low-mass microquasars, which could produce a significant 

Q.1 luminosity in the gamma-ray domain. 

Methods. We consider that a significant fraction of the relativistic particles in the jets of low-mass microquasars are 
protons and then we study the production of high-energy emission through proton synchrotron radiation and photopion 
production. Photopair production and leptonic processes are considered as well. We compute a number of specific models 

CZ3 , with different parameters to explore the possibilities of this scenario. 

Results. We find that important luminosities, in the range of 10 34 — 10 37 erg s _1 , can be achieved by proton synchrotron 
radiation in the Gamma-Ray Large Area Space Telescope (GLAST) energy range, and lower, but still significant 
luminosities at higher energies for some models. 

kJ Conclusions. We conclude that the proton microquasar model offers a very interesting alternative to account for the halo 

gamma-ray sources and presents a variety of predictions that might be tested in the near future by instruments like 
GLAST, the High-Energy Stereoscopic System II (HESS II), the Major Atmospheric Gamma-ray Imaging Cherenkov 
r"Si ' telescope II (MAGIC II), and neutrino telescopes like IceCube. 
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OO , 1. Introduction ral arms of the Galaxy. Imaging air Cherenkov telescopes 

,. „ s ,, „ . . . ,. have recently detected HMMQs (e.g., Cygnus X-l and LS 

• . Microquasars (MQs) are accreting X-ray bmanes that, pro- j +gl 303 b the Majop At heric Gamma-ray Imaging 

S C f re J latlvlstlC 1 P artlcles (e.g., Mirabel & Rodriguez cherenkov telescopej MAGIC, Albert et al. 2006, 2007; and 

1999). According to the type of donor star, microquasars Lg 5Q39 b the ffi gh-Energy Stereoscopic System, HESS, 

m n , Jf,n It ™^ g S ^1 l OW - masS microquasars Aharonian et aL 2005) There is some controversy around 

(HMMQs and LMMQs respectively). Soon after their ■ cbs- the nature of gome of thege objectgj especiall of LS r +61 

covery by Mirabel et al. (1992), it was proposed that MQs 3Q3 we refer the re&der tQ Romero ^ &{ (200?) for & dig _ 
might be sources ol gamma rays, as it is the case lor • 

some active galactic nuclei, in particular, fiat-spectrum ra- _ ' , , . „ ... . , , .„ . 

,. t nr T /-At t \ l • tr o 1 here is a population of variable unidentified sources 

dio quasars and HL Lacertae (HL Lac) obiects (Levmson & , . , , , , , ^ i , • r^ nnn , 

„, , r j in n^\ a r li .i j i j forming a halo around the Galactic center (Gremer 2001, 

Blandford 1996). A variety of models were then developed , T , ° , „„„„> T , . , ,, , ,, . v r 

c niinjA ii u j , u i Nolan et al. 2003). ft is expected that this group of sources 

lor HMMQs, mostly based on interactions between rela- ,, , / ,. , ., " , . .... °i i tu i i 

, . . ,. ,. , • ,i ■ n ,i i j, j \ j is old and presents a distribution similar to that of the glob- 

tivistic particles m the lets (both leptons and hadrons I and , , , , ... ,. -, , 

, , n !i r ,i i ; • • , ular clusters, but with no direct positional correspondence, 

matter and photon holds from the not, massive primary star n . ,, , . v% ., , ■, , j 

i -n onnr j r A • t> i r, Since they are extremely variable, it has been suggested 

(see Romero 2005, and references therein; Bosch-Ramon ,, , , . . i,, T inm /• u- j. \ n j 

\ , nnn „ t, , t, nnnrr i n ii • \ that they might be LMMQs (i.e., compact obiects) expelled 

et al. 2006; Bosch-Ramon 2007, and references therein). , ii.ii i., \ .■ i m ■ ™m 

m, . . . r ,, , , , / from globular clusters or the galactic plane (Gremer 2001, 

Ihe primary motivation tor these models was to explain „„„.! v 

the existence of a population of variable and unidentified _/* 

gamma-ray sources on the galactic plane (Romero 2001; , Dlfferent types of models have been proposed to ex- 
Kaufman Bernado et al. 2002; Bosch-Ramon et al. 2005). P^m gamma-ray emission m LMMQs, where the exter- 
rp, i ,ii , j • .i • • nal holds to the ict itself arc exiguous. Aharonian & 
I hese sources are young and mostly located in the inner spi- , ,„„„„% , . „ ,, . /,„„„n , 
6 J Atoyan (1998) and Atoyan & Aharonian (1999) developed 

Send offprint requests to: G. S. Vila time-dependent synchrotron self-Compton (SSC) models 
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model with emission reaching the hard X-rays. Grenier et 
al. (2005), on the other hand, calculated the interaction 
of a relativistic leptonic jet with the photon fields from the 
corona, the accretion disk and the dim star, concluding that 
external inverse Compton interactions cannot be invoked 
to explain the high gamma-ray luminosities (in the range 
10 34 — 10 36 erg s _1 ) inferred from the Energetic Gamma- 
Ray Experiment Telescope (EGRET) detections. 

Hadronic models based on inelastic proton-proton col- 
lisions (e.g., Romero et al. 2003) are not expected to work 
in the bare environment of old and cold stars as those of 
LMMQs. However, Levinson & Waxman (2001) proposed 
that photopion production could be important in MQs if 
strong synchrotron radiation fields are generated at the 
base of the jets. These authors developed a simple model to 
predict the neutrino emission of MQs, where the neutrinos 
are generated through the decay of charged pions created 
in proton-photon interactions. 

Here we explore in more detail a proton model for 
LMMQs, with emphasis on the electromagnetic output. We 
consider a system where the jets have a content in relativis- 
tic protons, which are accelerated by shocks in a compact 
region close to the accreting object. Electrons are acceler- 
ated as well in this region. We study models with different 
ratios of primary relativistic protons to electrons. The jet 
itself is considered as a mildly relativistic macroscopic fluid 
launched by some type of magneto-centrifugal mechanism, 
and hence it is initially assumed to be in equipartition with 
the magnetic field (as considered, e.g., by Bosch-Ramon 
et al. 2006). The jet power is taken as a fraction of the 
accretion power. In the case of advective accretion flows, 
the spectral energy distribution will be dominated by the 
non-thermal emission. Among the processes that cool the 
relativistic particles we consider synchrotron radiation for 
both electrons and protonqj, inverse Compton (IC) interac- 
tions with the locally produced fields, photopair and pho- 
topion production, and adiabatic losses. The results show 
that proton microquasars might be detected by forthcoming 
satellites like the Gamma-Ray Large Area Space Telescope 
(GLAST) and by ground-based air Cherenkov telescopes 
with a low-energy threshold. 

The structure of the paper is as follows. In Sect. 2 we 
outline the basics of the model and we introduce all the rele- 
vant radiative processes. In Section 3 we present the results 
of our calculations for a number of specific models. Section 
4 presents the effects of internal photon-photon absorption 
on the production spectra. Section 5 contains a discussion 
of the results. We close in Sect. 6 with a brief summary and 
the prospects of our model for future applications. 



disk-jet coupling hypothesis of Falcke & Biermann (1995), 
the kinetic power of each jet is assumed to be proportional 
to the accretion power, L ict = q iC tL&ccr, with q- ict — 0.1 
(Falcke & Biermann 1995; Kording et al. 2006). A fraction 
L TC i = Qrci^jct of this power is in the form of relativistic 
particles. We include both hadronic and leptonic content, 
L rc \ = Lp + L e . The way in which energy is divided between 
hadrons and leptons is unknown, but different scenarios can 
be taken into account by setting L p = aL e . In order to have 
a proton-dominated jet, we consider models with a > 1, as 
well as the case of equipartition between hadrons and lep- 
tons, a = 1. 
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Fig. 1. Top panel: Components of a low-mass microquasar. 
Bottom panel: Detail of the jet launching region with the 
relevant geometric parameters. 



2. Model 

2.1. Basic scenario 

The low-hard state of X-ray binaries is characterized by 
the presence of steady jets and a hot corona around the 
compact object (see Fig. [1]). We assume that the jets are 
launched perpendicularly to the orbital plane of the bi- 
nary, at a distance zq from the compact object and with 
an initial radius r^ = 0.1 zq. Afterwards the outflow ex- 
pands as a cone of radius r(z) = rp (z/zq). Following the 



1 The potential of proton synchrotron radiation in the context 
of AGNs has been explored by Aharonian (2002). 



The magnetic field in the jets decreases with the dis- 
tance z to the black holeQ. If the flow expands adiabatically, 
B(z) is given by 



B(z) 



Bo(^ 



(1) 



MQs with a weakly-magnetized neutron star are also possi- 
ble, e.g. Cir X-l and Sco X-l; see Migliari & Fender (2006) for 
a discussion. 
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The value of Bq = B (zq) can be estimated by equipartition 
between magnetic and kinetic energy densities at zq, 



B% 



L 



jot 



8tt 2irr£v ict 



(2) 



Here Vj Ct is the bulk velocity of the jets. We find Bq « 
2 x 10 7 G for a jet with parameters similar to those of the 
well-known LMMQ XTE J1118+480, see Table [TJ 



Table 1. Model parameters 



Parameter 



Value 



L accr : accretion power [erg s _1 ] 

q^t- disk-jet coupling constant 

g rc i: jet's content of relativistic particles 

a: hadron-to-lepton energy ratio 

zq: jet's launching point [cm] 

z max : extent of acceleration region 

Tjct: jet's bulk Lorentz factor 

-Bo: magnetic field at base of jet [G] 

a: particle injection spectral index 

r\: acceleration efficiency 

E m ' n : minimum particle energy 

9: viewing angle 



1.7 x 10 39(1) 

0.1 

0.1 

1 to 10 3 

10 8(2) 

5zo 
1.5« 

2 x 10 7 
1.5-2.2 
10 -4 -0.1 

2 - 100 mc 2 
30° 



(1) Typical value for the LMMQ XTE J1118+480. 

(2) 50.Rschw for XTE J1118+480 (8M Q black hole). 



i] = 10 , describing an efficient and a poor acceleration, 
respectively. In Fig. [21 we show the cooling rates at zq for 
different processes of energy loss, together with the accel- 
eration rate, for protons and leptons. Under the physical 
conditions of our model the main channel of energy loss for 
electrons is synchrotron radiation and we can approximate 
*ioss ~ ^synchr' f° r P r °t° ns both synchrotron and adiabatic 
losses are relevant, t^ oss ps t~ ynchT + t~ d . The synchrotron 
cooling rate for a particle of mass m, charge e, and energy 
£ in a region of magnetic energy density Ub is 



4 /m e y 
3 V m ) 



, 4 /"m e \ 3 cotUb E 

synchr = 3 y m 



2 ? ' 

m e c mc 



whereas 



.-1 _ 2 v jct 

ad 3 z ' 



The maximum energy is thus fixed by requiring 

t^ c (E,z)=t^ s (E,z), 

For electrons this yields 



1/2 
E : aax (z) w 1.3 x 10 10 ( 77— ) c\'. 



(6) 



(7) 



(8) 



(9) 



and for protons 



In our model, particle acceleration is constrained to take 
place in a compact region near the base of the jets, extend- 
ing up to a distance z max of a few times zq. We assume an 
injection function that is a power-law in the energy of the 
particles, 



•(E,z) = Q Q 



E- 



erg 



_1 cm 3 . 



(3) 



This type of spectral distribution arises naturally in astro- 
physical environments as a consequence of diffusive particle 
acceleration by shock waves. Both the cases of hard and soft 
particle injection were considered in our calculations. For 
these spectra we adopted a = 1.5 and a — 2.2, respec- 
tively. Also, we added an explicit z-dependence to reduce 
the injection rate from the distance to the compact object. 
The normalization constant Qo for each type of particle was 
obtained from L e „ as 



L e , P = I d"r 



E ,„ax (z) 



^^e,p -&e,p ^ce,p\-^e,p-! %)* 



(4) 



The maximum energy that a relativistic particle can at- 
tain depends on the acceleration mechanism and the differ- 
ent processes of energy loss. The acceleration rate t~^ c = 
E~ 1 dE/dt for a particle of energy E in a magnetic field B 
is given by 



t: 



r/ecB 
E ' 



(5) 



where rj < 1 is a parameter that characterizes the efficiency 
of the acceleration. Two cases were considered: rj = 0.1 and 



E™ ax (z) a 4.5 x 10 



ic 



Z() 



1/2 



eV. 



(10) 



Particles cool rapidly after leaving the acceleration region, 
as can be seen from Fig. [31 which shows the evolution of 
E max (z). The size of the acceleration region can impose 
constraints on the maximum energy; this occurs only in 
the case in which the particle giroradius T s (z) = E/eB{z) 
is greater than r- ie t(z). Nevertheless, in our model maximum 
energies are in all cases given by Eqs. ^ and (fTU)) . 

We have calculated the steady state particle dis- 
tributions N(E, z) in the "one-zone" approximation 
(Khangulyan et al. 2007). This approximation is valid if 
the losses are very strong in the acceleration region and 
diffusion can be neglected. Then the transport equation 
(Ginzburg & Syrovatskii 1964) can be written as 



d_ 

dE 



dE 
~dt 



N(E,z) 



loss 



N(E,z) 



t e 



= Q(E, 



(11) 



Here t csc is the particle escape time from the acceleration 
region, which we estimated as t csc sa -z max /vj Ct . This equa- 
tion includes energy losses, removal of particles, and injec- 
tion. It provides a one-zone description of the flow since, 
for each z, it neglects the contribution to N (E, z) of parti- 
cles coming from other regions in the jet. It remains a good 
approximation as long as the acceleration/emission region 
is compact (see Khangulyan et al. 2008 for a different ap- 
proach). This is better fulfilled in the case of leptons, since 
they cool almost completely near the base of the jet where 
the magnetic field is stronger. The exact analytical solution 
to Eq. (fTTj) can be found in Khangulyan et al. (2007): 
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Fig. 2. Acceleration and cooling rates at the base of the jet for protons (left panel) and electrons (right panel), calculated 
for representative values of the model parameters (E min — 100 mc 2 , a = 1 and a = 2.2). 
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Fig. 3. Evolution of the maximum (kinetic) energy of protons (left panel) and electrons (right panel) with the distance 
z to the compact object. Outside the acceleration region particles cool rapidly due to adiabatic and synchrotron losses. 



N{E,z) 



where 



dE 



dt 



-1 ,S m "( Z ) 
loss J E 



dE'Q(E',z) 

r(E,E') 



x exp 



r-E' 

t(E, E') = / dE" 
Je 



dE'' 



dt 



> (12) 



(13) 



of protons and electrons, proton-photon inelastic collisions 
(p7) and IC scattering. 

To calculate the synchrotron spectra we use the classical 
formulae, see for example Blumenthal & Gould (1970). In 
the jet co-moving reference frame, where particle distribu- 
tions are isotropic, the power radiated by a single particle 
of energy E, and pitch angle a is given by 



P.,.„ h , (E 7 , E, ., a) - «> f r di KtMQ 



E-y/E C 



loss 



(14) 



The resulting distributions have a power-law dependence on 
the energy of the particles. In the case of protons, N p (E p , z) 
mimics the injection function, N p (E p ,z) oc E~ a . The syn- 
chrotron cooling timescale is much shorter for electrons 
and they radiate a significant fraction of their energy be- 
fore escaping. This modifies the injection spectrum mak- 
ing it softer, and N e (E e ,z) can be well approximated as 

N e {E ei z)(X Ee 



where K 5 / 3 (() is a modified Besell function. The syn- 
chrotron power peaks sharply around the characteristic en- 
ergy E c , given by 



E r 



3heB(z)sma ( E 



(15) 



(cH-l) 



Denoting the volume of the emission region by V, the total 
luminosity is 



2.2. Radiative processes 

We consider three processes of interaction of relativistic 
particles with the fields in the jet: synchrotron radiation 



L. 



synchr 



{E n 



£L 



d 3 r / dSl Q sin a 
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'0) 



dEN(E,z)P synchr . (16) 



Electrons in the jet interact with the synchrotron ra- 



Proton-photon (jry) interaction has two main channels: 
photopair production 



diation fields through IC scattering. As leptons are highly .,, ,, , , , v 1 n/r ir -<• j-v. u *. 

,,..,.,, , . , , c . „,, . with a threshold energy 01 ~ 1 Mev lor the photon in 

relativistic they can lose an important traction ol their en- , -, , , c 

, ,,. . , ,, T ,, . TL.T- i • ;• the proton rest frame 

ergy in each collision and trie Klein- JNismna cross section , ,-, , , . , , 

OJ . , , . ,. , proceeds through two main branches 

must be used in the calculations. According to Blumcnthal 

& Gould (1970), the spectrum of photons scattered by an 

electron of energy E c = 7 c 77i c 2 in a target radiation field 

of density n(e, z) is 



p + "f^p + e + e + , (25) 

with a threshold energy of ~ 1 MeV for the photon in 
the proton rest frame, and photomeson production, which 

(26) 



Pic (E 7 ,E e ,e,z) 



3cttc n (e, z) 



F{q). 



(17) 



p + 7 — > p + an + b (n + + 7r~) 
and 

p + 7 — ► n + 7r + + air + b (ir + + 7r~ ) 



(27) 



with approximately the same cross section and a threshold 
energy of ~ 145 MeV. Here a and b are the pion multiplic- 
Here a T is the Thomson cross section and the function F (q) itieg charged pions decay in leptons and neut rinos, 



is given by 

F(q) = 2qlnq+(l + 2q) (1-q) 

1 









+^(1-9) 



(<?r c 



(i + r e g)' 



where 



r c = 

and 

q = 



4e7c 
m c c 2 



E~.: 



F C E C (1 — Ej/E e ) 



(18) 



(19) 



(20) 



whereas neutral pion decay yields photons 

7T° -> 2 7 . 



(28) 
(29) 

(30) 



The pr/ cooling rate for a proton of Lorentz factor 7 P in a 
photon field of density n (e) is 



-t ( \ 

t vl ,iVlp) = ^2 



_ra(e) 



2e 7p 



deV^^eO-K^i (e') e', 



(31) 



The limit T c -C 1 corresponds to Thomson scattering, but 

Eq. fiTj) is valid for all T e , even deep into the Klein-Nishina where i = e ± ,TT denotes the interaction channel, e' is the 
regime, as long as 7 e 3> 1. The allowed range of energies energy of the photon in the rest frame of the proton, and e t h 

is the photon threshold energy measured in the same frame; 
K V1 is the inelasticity of the process. The corresponding- 
collision rate is given by a similar expression: 



for the scattered photons is 

r. 



e< £ 7 < 



1 



-E r 



(21) 



Given the electron distribution TV (E e , z), the total IC spec- 
trum is calculated as 



i(7 P ) 



2 % J^ 



dc 



n(e) 



2«7p 



de'cTp 7i i (e')e'- 



(32) 



Lie (Ej) — E, 



d 3 ? 



B » a x (2) 



dE e N{E c ,z) 



To estimate the spectrum from the decay of 7r°, we followed 
Atoyan & Dermer (2003). In this approach, the cross section 
cr p7j 7t and the inelasticity -HT P7l -n are approximated by step 
functions 



dePic- (22) 



In our case, the scattered photons are those of the syn- 
chrotron radiation fields of both protons and electrons, 



r(eO 



-""P7, 7T \^ J 



n (e, z) = n synchr (e, z) . 



(23) 



340 /ibarn 200 MeV < e' < 500 MeV 

120/ibarn e'>500MeV, 
r 0.2 200 MeV < e' < 500 MeV 
I 0.6 e'>500MeV. 



We estimated n sync h r in the /oca? approximation of 
Ghisellini et al. (1985), 



"synchr v^' *J 



e sy (e,z)r(z) 



(24) 



These energy ranges correspond to the single-pion and 
multiple-pion production channels, respectively. 

The spectra are calculated in the (^-functional approxi- 
mation in the energy of pions and photons. In the single- 
pion production channel, each tt° is produced with an en- 
ergy rj K\E p and this energy is equally distributed among 



where e synchr (e, z) is the synchrotron power per unit the products of its decay, yielding 



volume per unit energy. 



E-, 



1 



\KiEp = 0.1£ p 



(33) 
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In the multiple-pion production channel almost all the en- 
ergy lost by the proton is equally distributed among three 
leading pions 7r°, 7r + , and ir~ . The mean energy of each 
pion is « KiE v fi and, therefore, the energy of the photons 
results the same as in the single-pion production channel. 
If pi and p2 = 1 — pi are the probabilities of the pj 
collision taking place through the single-pion and multiple- 
pion channel, respectively, and £ pn ~ 0.5 is the probability 
of conversion of a proton into a neutron per interaction, the 
mean number of neutral pions created per collision is 



n„o =pi(l- C pn ) +P2- 

We can define a mean inelasticity 

K P1 ^ (7 P ) = t~ (7p)w P7iTr (7 P ) , 



(34) 



(35) 



then the probabilities pi.2 can be calculated from the rela- 
tion 

K P1 ^{E p )=p 1 K 1 + {l-p 1 )K 2 . (36) 

The emissivity of ir° in the ^-functional approximation is 



Q-n (Ejr) — / dEpNp (E p )u> pl ^ (Ep) 

xn»o (E p ) 5 (E„ - 0.2E P ) 
= 5N P (5£V) lo p1 ,k (5E W ) n^o (5E n ) . 



(37) 



Taking into account that each 7r° gives two photons, the 
photon emissivity results 

Q 7 [EL,) = 2 J AE n Q n a (E„) S (£ 7 - 0.5£ w ) 

= 20N P (1QEJ uj pi ^(\0E 1 ) n„o (10 E^) . (38) 

The emissivity of secondary pairs can be estimated in the 
same way. Both in the single-pion and in the multiple-pion 
channel each charged pion has an energy ~ 0.2E p . This 
energy is equally distributed among the products of meson 
decay, hence the energy of each electron/positron is, on 
average, E e ~ 0.05E p . The mean number of charged pions 
created in a p^ collision is 



n*± = £,pnPl + 2P2, 



(39) 



and since only one lepton is produced in each decay, the 
emissivity of pairs is given by 



Q e ± (E e ±) = 20N p (20E e ±)u> I ry,i r (20E e ±) 

xn I± (20E e ±). (40) 

The second channel of proton-photon interaction is pair 
production. This process has been studied, for example, 
by Chodorowski et al. (1992), Miicke et al. (2000) and 
Mastichiadis et al. (2005). The emissivity of pairs in the 
5-functional approximation is calculated as for the previ- 
ous process. Now the inelasticity can be approximated by 
its value at the threshold, K p7<e ± = 2m e /m p . Therefore, 



dEpNp (E p ) Wp7ie ± (E p ) S[E e ±- -^E p 

^E e A , (41) 



m p / mp \ I m ri 

2—Np —E e ± uj pl . e ± 
nip \ m P I \ Top 



An useful parametrization of the cross section a pl e ± can 
be found in Maximon (1968). 

All the expressions above are valid in the jet co- 
moving reference frame, where the particle distributions 
are isotropic. The spectral energy distributions in the rest 
frame of the observer can be obtained by applying the ap- 
propriate boost. In what follows, we denote by primed and 
non-primed quantities those measured in the jet co-moving 
frame and in the observer frame, respectively. Photon en- 
ergies in both frames are related by the Doppler factor D, 



E 1 = DE' 7 . 



(42) 



The factor D depends on the jet's bulk Lorentz factor Tj C t 
and on the viewing angle 9. For an approaching jet, 



D 



1 



Fjet (l-/3jetCOS( 



(43) 



Finally, the luminosity in the observer frame is given by 
(e.g., Lind & Blandford 1985) 



i 7 (EL,) = D 2 L' (EL 



7 V 7/ 



(44) 



None of our results depend strongly on the value of the 
viewing angle as long as is it is not very close to zero, and 
we fixed 8 = 30° for calculation purposes. The values of the 
different parameters of the model used in the calculations 
are shown in Table Q] 

3. Results 

Figures 0] and [5] show the spectral energy distributions 
(SEDs) calculated for different values of the model param- 
eters. In Fig. Q] the spectral index of the injection function 
has the canonical value a = 2.2, whereas in Fig. [5] we con- 
sidered a harder injection with a = 1.5. For each case, we 
show four representative spectra obtained by varying the 
value of the remaining parameters within a physically rea- 
sonable range. The main contributions to the SEDs are al- 
ways due to synchrotron radiation of leptons and protons. 
The relative importance of each contribution depends on 
the value of the ratio a — L p /L e . The proton synchrotron 
spectra are hardly affected by changes in a, peaking at 



~ io J °- J6 erg s _1 for E 1 ~ 10 8_J eV; lepton synchrotron 
luminosities range from ~ 10 34 erg s _1 for a completely 
proton-dominated jet (a — 10 3 ) to ~ 10 37 erg s _1 in the 
case of equipartition (a = 1). In fact, lepton energy losses 
are so strong that they radiate almost all their available 
energy budget in each case. 

The efficiency of the acceleration r\ fixes the maximum 
particle energy and therefore has a direct effect on the high- 
energy cutoff of the SEDs. Synchrotron spectra extend up 
to E 1 ~ 10 12 eV in the case of protons and up to E~, ~ 
10 9 eV in the case of leptons for r\ = 0.1; when a poor 
acceleration efficiency r\ = 10~ 4 is considered, only energies 
of about three orders of magnitude smaller are reached. 
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Fig. 4. Spectral energy distributions obtained for different values of the parameters in the case of a particle injection 
index a = 2.2. 



In the same way, modifying the minimum particle energy 
accordingly changes the low-energy cutoff of the spectra. 

The IC contribution is always negligible and, for a > 1, 
IC luminosities are below ~ 10 30 erg s _1 as the leptonic 
content of the jet is reduced. Photopion production (^7^0) 
yields a very hard energy tail to the SEDs, peaking at en- 
ergies E 1 ~ 10 14-15 eV. The highest value of the p^^ lu- 
minosity is sensible to a, since the target photons for p-f 
collisions are provided mostly by the lepton synchrotron ra- 
diation field. Luminosities as large as 10 35-36 erg s _1 can be 
reached at very high energies. Synchrotron radiation from 
e produced in pj e ±^± can be important in those models 
with relatively low barion content, as it can be seen from 
the different SEDs in Figs. H] and In general, these con- 
tributions have luminosities not too different from those 
obtained from pj„o but they cover a lower energy range. 



4. Internal photon-photon attenuation 

Internal photon-photon absorption might be important in 
some of the models we have calculated. In particular, TeV 
gamma-rays can be efficiently absorbed by IR synchrotron 
photons. This might initiate an electromagnetic cascade in 



the emission region of those cases where there are signifi- 
cant photon fields at low energies as could happen in blazars 
(e.g., Blandford & Levinson 1995). In the context of micro- 
quasars, however, a strong magnetic field can quench a pure 
IC cascade. This is the result of synchrotron losses for the 
first few generations of pairs, which then cannot produce 
photons energetic enough to sustain the electromagnetic 
avalanche (Khangulyan et al. 2008). The photon absorp- 
tion, nonetheless, might change the gamma-ray spectrum 
depending on the specific model. This, in turn, can yield a 
variety of slopes in the spectra to be observed by GLAST 
or other instruments sensitive in the MeV-up-to-TeV en- 
ergy range (see, e.g., Aharonian et al. 2008 for a recent 
discussion in the context of extragalactic jets). 

The opacity of a photon field of density n(e, z) given 
by Eq. (J24j) to the propagation of a gamma-ray photon of 
energy E n is (Gould & Schreder 1966): 



r(E 7 ) = - 



l •'eth 



(1 — it) <r(Ey, e, u) 



xn(e,z) du de dl, (45) 



G.E. Romero and G.S. Vila: The proton microquasar 



£ min =2mc 2 - o=1 .5 - n=0.1 - fl=1 



£ min =2mc 2 -a=1.5-T l =0.1 -a=10 3 



_— 35 
[/) 

£, 34 

_1 

g> 33 

_l 

32 



synchrp 


— • — synchr photopion ■ 


synchr e- ^_^_ 


_^ —* — synchr photopair - 


IC ^^ 


\ 




\ ' \ /\ 


- Pv ^^ 


v , •* 1 •■ 




/\ y i •• ■• 


1 / 




1 / 
1 ' 




: / 

■ / // 


A 

A 




i|; 


1 f ' 


A 


/ •', 


1 


7 ' ' 


A 




. ■ 


/ '/ 


t 


/ • ! 


: 




A, 1 L .1 . 1 . 1 





10 12 14 16 



37 

36 

,=— 35 
w 
O) 

& 34 

_i 

g> 33 

_l 

32 



- synchr photopion 
-synchr photopair 



■synchrp 
- synchr e- 

IC 

Pi 




10 12 14 16 18 



Log E (eV) 



Log E (eV) 



£ rain =100mc 2 -a=1.5-T l =0.1 -a=100 



E^IOOmc 2 - o=1 .5 - rpIO" 4 - a=100 



.r- 35 

ai 

& 34 
_T 

o> 33 



- synchr photopion 
-synchr photopair 




Log E (eV) 



j*" 33 



-•— synchr photopion 
-• — synchr photopair 




- synchr p 
- synchr e- 
IC 
■Pi 



\A 




1 

1 
1 


A 

A: 



Log E (eV) 



Fig. 5. Spectral energy distributions obtained for different values of the parameters in the case of a particle injection 
index a — 1.5. 



where u = cosd, ■d is the angle between the momenta of 
the colliding photons, I is the photon path, and the cross 
section for the interaction is given by (e.g., Levinson 2006): 



,(25 71 e, 0) = -cr T (l - 
3 - /3 4 ) In 



/3 2 ) 

'1 + 

1-/3 



- 2/3 (2 - p 2 ) 



(46) 



In the latter expression, /3 is the speed of the elec- 
tron/positron in the center of momentum frame. It is re- 
lated to the energy of the incident and target photons by: 



(1-/3 2 ) 



(1 -u)E^e 



0</3< 1. 



(47) 



1 



The threshold energy e t h is defined by the condition /3 
with § = 0. 

We have calculated the attenuation factor exp(— r) for 
the different models presented in the previous section. Some 
results are shown in Fig. [6] 

We notice that in all cases with significant leptonic syn- 
chrotron emission, the radiation is completely suppressed 
above 10 GeV. For proton-dominated models, which are 
characterized by a prominent proton synchrotron peak, the 
attenuation is quite moderate. 
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Fig. 6. Attenuation curves for different models of a proton 
microquasar. In all cases the acceleration efficiency is r\ = 
0.1. 



In Fig. [7J we present two spectral energy distributions 
modified by the effects of the photon-photon absorption. In 
the left panel, we show a case with a = 1 (strong atten- 
uation). The internal opacity to gamma-ray propagation 
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results in a significant softening of the spectrum between 
~ 10 MeV and ~ 10 GeV, with a cutoff beyond the latter 
energy. This is consistent with the type of very soft spectra 
observed by EGRET in some variable halo sources (e.g., 
Thompson et al. 2000, Grenier 2001). In the right panel of 
the same figure, we show the SED for a proton-dominated 
microquasar, which is basically not modified by absorption 
effects. 

In all cases the emission at TeV energies is either sup- 
pressed or too weak, and hence no detection with current 
Cherenkov telescope arrays should be expected for these 
objects unless very high acceleration efficiencies could be 
achieved in proton-dominated casetQ Neutrino propagation 
is not affected by the effects discussed in this section. 

5. Discussion 

From the results presented in the previous sections is 
clear that proton low-mass microquasars can be significant 
gamma-ray sources in the MeV-GeV range, and perhaps in 
some cases even at TeV energies. According to the ratio 
a of primary protons to leptons, we can obtain different 
types of low-energy counterparts. Such counterparts range 
from sources with radio luminosities of 10 29 — 10 30 erg s^ 1 
and strong X-ray emission (e.g., the model shown in the 
left upper panel of Fig. [5]), up to sources with weak lumi- 
nosities dominated by proton synchrotron emission between 
100 MeV and 100 GeV (e.g., the model shown in the right 
upper panel of Fig. [5]). Since the cutoff of the synchrotron 
peak always occurs at energies below 1 TeV, and the high- 
energy emission due to photopion production is beyond 10 
TeV in all models, detection with ground-based Cherenkov 
telescopes with a threshold of around 1 TeV would be dif- 
ficult. In cases with strong leptonic synchrotron emission, 
the TeV radiation is directly suppressed by internal ab- 
sorption. The telescope GLAST, on the contrary, with its 
energy window of 100 MeV — 300 GeV is especially suit- 
able for the detection of these objects. Actually, we suggest 
that many of the unidentified EGRET sources detected off 
the galactic plane around the Galactic center, might be 
proton microquasars. Variability, one of the outstanding 
properties of these sources (Nolan et al. 2003), can be eas- 
ily introduced in our model through a variable accretion 
rate, precessing jets, or internal shocks (in this latter case, 
there would be very rapid variability, superposed to longer 
variations with timescales from hours to days). The Astro- 
rivelatore Gamma a Immagini LEggero (AGILE) satellite 
might provide a more accurate location for some of these 
sources in the immediate future. Low-threshold, ground- 
based gamma-ray Cherenkov telescope arrays like HESS 
II and MAGIC II could also detect hadronic LMMQs at 
E > 100 GeV. In particular, the tail of the proton syn- 
chrotron peak from models with large values of a might 
be detectable, displaying a soft spectra. So, observations of 
LMMQs with Cherenkov telescopes can be useful to con- 
strain the parameter a. 

An additional prediction of our model is the produc- 
tion of high-energy (E > 1 TeV) neutrinos, with luminosi- 
ties in the range 10 33 — 10 35 erg s _1 . In the cases of the 
highest luminosities, if the MQ is located not too far away 
(say, around 2 kpc as it is the case of XTE J1118+480) 



3 A high acceleration efficiency would move the high-energy 
cutoff of the synchrotron peak into the TeV regime. 



the expected flux could be similar to those estimated for 
HMMQs (e.g., Levinson & Waxman 2001; Romero et al. 
2003; Christiansen et al. 2006; Aharonian et al. 2006). We 
notice that heavily absorbed sources are precisely those 
with the highest neutrino luminosities, so a correlation be- 
tween gamma-ray and neutrino fluxes should not necessar- 
ily be expected. Finally, we mention that some features re- 
lated to electron/positron annihilation could be observable 
in cases with strong absorption (e.g., Romero 1996; Punsly 
et al. 2000; Levinson 2006). 

6. Conclusions 

We have developed a simple model for proton low-mass mi- 
croquasars that predicts significant gamma-ray emission in 
the MeV and GeV domain for these objects. Observations 
with GLAST will help to establish some fundamental pa- 
rameters such as the ratio of relativistic protons to elec- 
trons in the sources and the effects of internal absorp- 
tion. Proton low-mass microquasars also produce very high- 
energy gamma-rays (well within the TeV range) and neu- 
trinos, that in some cases might be detectable by km 3 - 
detectors like IceCube or Nemo. Proton microquasars might 
also present very weak low energy counterparts, and hence 
they are suitable candidates to explain the unidentified 
gamma-ray sources detected by EGRET in the galactic 
halo. 

Future developments of the model should include a more 
sophisticated treatment of the acceleration region and time- 
dependent calculations. These topics will be included in a 
forthcoming communication. 
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